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We relate the intermittent fluctuations of velocity gradients in turbulence to a whole range of local 
dissipation scales generalizing the picture of a single mean dissipation length. The statistical distri- 
bution of these local dissipation scales as a function of Reynolds number is determined in numerical 
simulations of forced homogeneous isotropic turbulence with a spectral resolution never applied be- 
fore which exceeds the standard one by at least a factor of eight. The core of the scale distribution 
f~ — ■ agrees well with a theoretical prediction. Increasing Reynolds number causes the generation of ever 

c 2 ^ ' finer local dissipation scales. This is in line with a less steep decay of the large-wavenumber energy 

\ spectra in the dissipation range. The energy spectrum for the highest accessible Taylor microscale 

04 ■ Reynolds number Rx = 107 does not show a bottleneck. 
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I. INTRODUCTION 



Turbulence is characterized by large fluctuations of velocity gradients which appear preferentially at the smallest 
scales of the flow. The arnplitudes of these events exceed the mean values by orders of magnitude which is known as 
small-scale intermittency [l], [2] ■ It is also believed that these intensive fine-scale fluctuations are intimately connected 
with the nonlinear cascade-like transfer of kinetic energy through the hierarchy of eddy structures that fill the fluid on 
^ \ larger scales [1, 0, H S 0] ■ ^ better understanding of fiuid turbulence as a whole requires thus a detailed resolution 
M-H of the intermittent dynamics at the small-scale end of the inertial range. In other words, it is necessary to determine 
c/2 how deep into the beginning dissipation range the roughest filaments from the inertial range can sweep and how they 
. ^ affect the spectral decay of turbulent fluctuations with increasing Reynolds number. This can also be important for 
the mixing in reacting and non-reacting flows at high Schmidt number Q where a significant fraction of stirring 
of the concentration field takes place in the dissipation range of the flow. Such a study requires that the steepest 
gradients and their statistics are well resolved. Although significant progress in measurement techniques has been 
made [10], the finest structures remain still spatially unresolved in experiments. 

In this Letter, we present numerical simulations of forced homogeneous isotropic turbulence that unravel the 
dynamics in exactly the region where inertial and dissipation ranges match, both in physical and Fourier space. 
In order to study the fine-scale structure and their statistics, a grid resolution is applied which exceeds that of 
standard simulations by an order of magnitude. In particular, the following questions will be answered. What are 
the finest spatial scales across which large-amplitude gradient events evolve? How are these scales distributed as 
■ a function of the Reynolds number? Is consequently the large-wavenumber behaviour of the energy spectra in the 
' dissipation range Reynolds number-dependent? 
I . The large-wavenumber behaviour of the turbulent fluctuations at the smallest scales is a long-standing problem. 
' Kolmogorov postulated 65 years ago a universal form of the decay of the energy spectrum that goes deep down into 
, the dissipation range ^V^. Since then several analytical attempts have been made to determine the spectrum. The 
IL^ ■ works left unspecified constants [H, [H, [13, [ip or considered an infinitely extended range of excited scales 0, [T^l ■ 
Numerically, the time advancement in the dissipation range is very challenging since a significant fraction of the 
computational ressources has to be spent for the resolution of the tiny amplitudes [H, [l^, [Ifl] . 



H , The classical theory of turbulence predicts a mean scale at which the turbulent cascade ends and the fiow viscosity 
. . . 1 starts to dominate. This scale is known as the Kolmogorov length [T]| 

m = j^, (1) 

where ly is the kinematic viscosity of the fluid. It is derived by a dimensional estimate and does not capture for the 
observed intermittency of the velocity gradients. We see that the energy dissipation rate, that probes the velocity 
gradient magnitude, enters the definition of rjK as a mean, (e). On the one hand, there are very intense gradients 
in the form of thin and stretched vortex tubes that seem to have diameters around r]K, or even less. On the other 
hand, ambient regions will exist with typical spatial variations larger than rji^. A whole range of local dissipation 
scales around the classical Kolmogorov length follows from the picture. This idea was put forward first within the 
multifractal formalism [3, 4]. However, a direct analysis of the scales from data records remained very difficult, simply 
because these structures were not resolved. 
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Yakhot derived an estimate that connects a dissipation scale rj with a velocity increment u.,, = \u{x + t]) — u{x)\ 
across this scale from the equations of motion (see [2l[ and references therein) . It reads 

r]u,j Ki u . (2) 

Technically, eq. ^ is derived from a local kinetic energy balance by a so-called point-splitting procedure [l^, [H, H^l ■ 
Equation ^ tells us that is a field fluctuating in space and time and suggests an implicit way to determine the local 
dissipation scales 77 in numerical simulations, by velocity increments over Kolmogorov and sub-Kolmogorov distances. 
Relation ([2]) can be also obtained by equating convective time scale, Tj/uy, and dissipation time scale, if' /v, a step 
which is at the core of the multifractal approach of Paladin and Vulpiani [j| . 

The probability density function (PDF) of the local dissipation scales will be denoted as Q{r]) in the following. In 
fig. [H we show an instantaneous snapshot of two isolevel sets of the field 77 where the isosurfaces in red are nested in 
the transparent gray ones. The figure clearly underlines the fiuctuating character of the field 77. 

Once one accepts the concept of local dissipation scales, the question on the smallest local dissipation scale arises. 
This scale can be determined by matching the inertial and dissipation range dynamics in the equations for the 
longitudinal increment moments of order 2n which are defined as 

S2n{r) = (u^) = {Hx + r) - «(x)|2") (3) 

where u is the turbulent velocity field projected onto r and r — |r| [7|,|23|. The matching scale becomes order-dependent 
and is given by 

r]2n = LR^'-^^-+'^\ (4) 

Here, L is the integral scale and R = glL/v the corresponding large-scale Reynolds number with ul = {uj^y^^- The 
exponents ^„ determine the scaling behaviour in the inertial range of turbulence, i.e. Sn{r) ~ r^" for rjK < r < L. In 
the limit of n — + 00 we expect C2n ~ Czn+i and thus 

Vmin = LR^^ . (5) 

One obtains a steeper decrease with respect to the Reynolds number as for the Kolmogorov scale, r/K — LR^^/^. The 
same relation ([5]) follows in the multifractal approach by Paladin and Vulpiani [3| for the roughest increments with 
Holder exponents h — 0. 



II. NUMERICAL MODEL 



The three-dimensional Navier-Stokes equations for an incompressible flow are solved. The direct numerical simula- 
tions are based on the pseudospectral method with fast Fourier transformations and a 2/3 de-aliasing. The simulation 
domain is a periodic cube with a volume ofV — {2n)^. The velocity field is kept in a statistically stationary state by a 
large scale volume forcing that is added to the r.h.s. of the Navier-Stokes equations. It injects kinetic energy into the 
flow at a fixed rate, em, and thus prescribes the mean energy dissipation rate. In the case of statistical stationarity 
the kinetic energy balance reduces to = (e). More details are found in 0]. The turbulence is homogeneous and 
locally isotropic. Some parameters of the present simulations are listed in table 1. The standard spectral resolution 
is determined by the criterion kmaxVK > 1-5 [25]. It is seen from the table that our applied resolution is at least by a 
factor of eight better than the standard case. 



III. DISTRIBUTION OF LOCAL DISSIPATION SCALES 



The derivation of Q{r]) starts with relation ([2]) which connects velocities and scales at the small-scale end of the 
inertial range ^2l\. Q{rf) can then be calculated from the PDF of the longitudinal velocity increments across distances 
•q which is conditioned on 7/^77/1/ = c (see eq. (H))). c is a constant 0(1) and it follows 

The PDF of the velocity increments is obtained from a Mellin transform [2^. In Ref. [l^], Mellin transforms were 
used for the first time in turbulence to construct PDFs from increment moments by 
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TABLE I: Parameters of the direct numerical simulations R\ — ^ 1^/ {{e)v)cr\ is the Taylor-microscale Reynolds number, and 
R = glL/v the large scale Reynolds number with the integral scale L and (t| — {{ul)^) (see eq. Q). The spectral resolution 
is indicated by kmaxijK where kmax = V2N/3. The mean energy dissipation rate (e) is 0.1 for all cases. 



Run 


N 


V 


L 


Rx 


R 




1 


512 


1/30 


1.02 


10 


12 


33.6 


2 


1024 


1/75 


0.92 


24 


32 


33.6 


3 


1024 


1/200 


0.76 


42 


74 


15.9 


4 


1024 


1/400 


0.69 


65 


143 


9.6 


5 


2048 


1/400 


0.69 


64 


140 


19.2 



6 2048 1/1000 0.66 107 347 9.7 




FIG. 1: (color online) Spatial distribution of local dissipation scales. An instantaneous snapshot of a simulation with 1024^ grid 
points is shown (see Run 4 in table 1). Isosurfaces are plotted for two levels, = Atik/S (shaded red) and 5riK /3 (transparent 
gray) . The isosurfaces for the smaller i] isolevel are nested within those for the larger one and are significantly sparser distributed 
in space. The figure underlines the fluctuating character of the local dissipation scale field. 



Three points are necessary in order to make progress. Firstly, prefactors in the scaling laws for the increment moments 
are fixed by recognizing their Gaussian statistics at the largest scale of the flow, L. Secondly, the inertial cascade 
range scaling law for the increment moments is still valid at the small-scale end of the inertial range such that we can 
write 

{uf) = {2p-l)Uaf{iy'\ (8) 

Thirdly, the unknown anomalous scaling exponents C,2p are approximated well with the polynomial C,2p = Sap — Abp^ 
for the lowest orders, p < 10 [2ll[23|. One finds a = (1 + 96)/3 = 0.383 and b — 0.0166 in order to satisfy the exact 
relation ^3 = 1. The Mellin transformation integral is then evaluated by a saddle point approximation and results for 
a given large-scale Reynolds number R to the following indefinite integral 

1 /■+°° 

Q{ri) — = / dxexp 

nriy/blog{L/r]) J -00 



(log( 



Ab\og{L/ri) 



(9) 



which can be analyzed by numerical quadrature. The distribution is supported for scales < rj < L only. 
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FIG. 2: Comparison of numerical and theoretical results from (21<] (dashed lines) for the probability density function of the 
local dissipation scale field at a given Reynolds number, Q{ri). Data are for Runs 1 (circles), 4 (asterisks) and 6 (squares) 
as given in Tab. 1. The values for c are 2.6 (Runl), 4 (Run4), and 4 (Run6) respectively (see Eq. (O). All other parameters 
in DNS and ((9| are identical. For a better visibility, the data for Runs 4 and 6 are shifted upwards by one and two orders of 
magnitude, respectively. Inset: Comparison of Q{r]) for the standard grid resolution with kmaxTiK = 1-2 and A*' = 128 (squares) 
and the present very-high-resolution case (solid line). Data are for run 4. 



The calculation of the PDF Q{ri) from the simulation data works as follows. A scale £ is fixed as an integer multiple 
of the grid spacing A, i.e. i — nA. The longitudinal velocity increments with respect to i in all three directions are 
determined at each grid site. If luijv k, 1, the grid site contributes to Q{t). The resulting distributions are shown 
in fig. O The inset underlines the necessity for the huge spectral resolution applied here. We compare a standard 
resolution case with the present one (which is eight times larger) while leaving all other parameters the same. The 
whole left tail of the local dissipation scale distribution cannot be resolved in the standard resolution case. 

The main picture of the figure compares our data for Runs 1,4 and 6 (see table 1) with the theoretical prediction 
from ([9|) . The distributions are rescaled by the scale 770 that arises from eq. jl]) when inserting = 2ap — 46p^ to 



since a = 0.383. We see that the distributions coincide quite well in the core and for most of the right tail of the PDF 
with the analytical shape. As stated in the caption, the parameter c becomes independent of the Reynolds number 
R which supports the use of the scale 770 for the rescaling. The smallest Reynolds number case (Run 1) however 
deviates. The reason might be the Gaussian statistics of the velocity gradients which changes to non-Gaussian for 
R\ > 10 — 15 (see Ref. [7| for a detailed investigation on this subject). 

Furthermore, we observe that the distributions start to deviate at the largest scales from the analytical shape for the 
two larger Reynolds numbers. This is attributed to the periodic boundary conditions in the simulations which affect 
the velocity increments taken over large distances. We have checked this by going up to velocity increments across the 
whole box length. Nevertheless, the slope of the algebraic decay for scales r\ > rjK does not vary significantly with the 
Reynolds number, neither for the data nor for the theory. This part of the PDF which corresponds with increments 
over larger distances remains almost insensitive to an increase of Reynolds number. Note that all data collapse there. 
They have been shifted in the figure for a better visibility. Stronger deviations arise in the left tail, i.e. for the finest 
scales. The analytical prediction is limited here. A next step would be to include higher-order corrections to the 



The numerical data for the left tail show a slight Reynolds number dependence. The value Qirj/rjo — 0.6) grows 
from approximately 10~^ at R\ = 65 to 10~^ at Rx = 107. The scales for the smallest Reynolds number go barely 
below 770- This indicates an increasing probability of very fine sub-Kolmogorov scales to appear. It is in line with 



IV. COMPARISON WITH NUMERICAL RESULTS 



rjo ~ LR-— ^ LR 
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FIG. 3: (color online) Compensated energy spectrum for Run 6. The red lines are 183 individual energy spectra saved during 
the time integration which took 4 large scale eddy turnover times, T^ddy = {u^)/L. The black line is the mean of all these 
spectra which is used for the subsequent analysis. The dashed line is drawn at 2. The inset magnifies the range of wavenumbers 
that should show a bottleneck effect. 




FIG. 4: Rescaled dissipation spectra for different Reynolds numbers. Wavenumbers are normalized by kp for which k^E{k) 
becomes maximal. The dotted line box indicates the data window from [2^ which is shown in the inset. The ratios kprjK 
were found to 0.167, 0.179, and 0.160 for Rx = 42,64, and 107, respectively. The numerical precision of the DNS at such 
small amplitudes was tested by a verification of the viscous decay law of single Fourier modes without nonlinear advection, 

Uk(i) = -!/fc^Uk(t). 



increasing small-scale intermittency of the velocity gradients which has been studied in [7| . The small-scale end of the 
support in each of the PDFs is taken as the scale "q^^^^ ■ As discussed in the introduction, the theoretical model 0,|2l| 
predicts rjmin = LR~^ (see eq. ([5])). Table 2 compares this estimate with our findings. We see that rj^'j^^ decreases 
more slowly than the theoretical prediction. The value of rjmin reached O.GriK for the largest Reynolds numbers. To 
conclude, ever finer sub-Kolmogorov scales are excited for increasing R, but not as pronounced as predicted by theory. 

The first moment of the distributions Qiij) gives a mean dissipation scale which is always larger than rjx- We find 
(rj) = 2.5, 6.2 and 8.1 for R\ = 10,65 and 107, respectively. For the range of accessible Reynolds numbers we thus 
observe an increase of the mean towards a scale ~ lOriK which is the generally expected crossover scale from the 
inertial to the viscous subrange [28j . 



FAR-DISSIPATION RANGE ENERGY SPECTRA 



The central question that was addressed in [T^, [T^ [13, is if the increasing small-scale intermittency in physical 
space manifests as well for the spectra in the crossover to the dissipation range or even in the far-dissipation range. 
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TABLE II: Ratio of the finest local dissipation scale -qmin to the box size ~ 2n following from the DNS data and from the 
estimate in i21| . respectively. is taken in order to show that /L^ is resolved, i.e. larger than l/N^. Furthermore, it 

provides an independent scale being the same in all runs. 



Rx 


10 


24 


42 


64 


107 


DNS 
' Irmn 

T}min — LFt 


1 

51 
1 

74 


1 

131 
1 

218 


1 

407 

1 
613 


1 

576 
1 

1295 


1 

1194 
1 

3232 



Figure [3] shows the instantaneous compensated energy spectra for the run at the largest Reynolds number. It provides 
information about the variation of the spectra in the dissipation range and how this adds up to a mean spectrum 
that will be used for the following analysis. Many numerical studies of the energy spectra have been focussed on 
the bottleneck phenomenon (see e.g. [30|). It is thought to result from a depletion of the nonlinear Fourier mode 
interaction at higher wavenumbers such that non-local mode couplings become more dominant. Wc did not observe 
a bottleneck for our spectra (see inset of fig. [3]). The reason can be the moderate Reynolds numbers which are 
accessed here. It can however also be that the significantly larger range of resolved Fourier modes in the viscous range 
diminishes this effect. 

Since the moderate Reynolds numbers prevent the quasi-algebraic sca ling analysis from Q , we study the spectral 
decay with respect to measures that follow from the spectra themselves |29j . The beginning of the dissipation range 
can be set at fc = fcj, where k'^E{k) has a maximum. In fig. [4l the dissipation spectra for the three largest R are 
rescaled by the corresponding kp. The dotted line indicates the data window from [2^. While the energy spectra 
collapse quite well within the dotted box, differences arise for the largest resolved wavenumbers. This indicates a 
deviation from the universality postulate made in the classical theory for the large-wavenumber decay of the energy 



spectra 11|. We wish to stress that our Reynolds number might still be too small for a firmer conclusion. 

The results from fig. 2] suggest a closer analysis of the Reynolds number dependence of the decay of the spectra in 
the far-dissipation range. The following form of the spectral decay for k — kr/K > 1 is used [l^ [ij, [H, 0, 

^(^) = ^5/^1,^)1/4 = PCk) = exp(-/3fc) (11) 

where a and P are Reynolds-number-dependent dimensionless constants. The exponential term can be motivated 
from Stokes eigenfunctions in the viscosity-dominated regime of the Navier-Stokes dynamics [1^. The power law 
term k" was calculated by the direct interaction approximation (DIA) to a = 3 [Tsj . Numerical experiments at low 
Reynolds numbers demonstrated later that a can exceed 3 due to small-scale intermittency which is not contained in 
the theory [l^. Equation pT|) is transformed into the following local slope form 

'-}2im^a-Pk (12) 
dlog(A:) 

which allows for direct determination of the two constants a and /3 by a least square fit. The upper picture of fig. [5] 
shows the local slope of the large-wavenumber spectral decay. The two lower pictures in the same figure list the 
results for both coefficients, aioc and /3;oc- We have performed therefore local least square fits to (IT2|) in the interval 
[fc — 1/2, fc -t- 1/2]. It should be stressed once more that the present data allow a systematic study of the dissipation 
range decay over an order of magnitude of Taylor microscale Reynolds numbers. The results are consistent with the 
findings from Refs. [ll[li,[23]. The exponent a, which determines the nature of the singularity in the Euler case, 
remains smaller zero for R\ > 24 which results in a convex shape of the dissipation range spectrum. Interestingly, this 
is the range of R\ for which velocity gradients are clearly non-Gaussian Q. The overall magnitude of the exponent 
P decreases for growing Reynolds number. A saturation of this decrease can be detected for the two largest Reynolds 
number runs. This would be consistent with a saturation to a constant magnitude which was predicted by Kraichnan 
p~3i |. Neither for a nor for /? a systematic behaviour with respect to wavenumber k is observed. The limited range 
of sub-Kolmogorov scales that can be resolved is one reason. Moreover, it should be stressed that the local fits of 
the exponential decay require a very high numerical accuracy as has been demonstrated in Pauls et al. [3l|. The 
double precision floating accuracy which we can provide only was not sufficient to apply more sophisticated local fit 
procedures to pT|) [sil |32| . Finally, it was also checked that the rescaling with 770 instead of rjK does not alter the fit 
results significantly. 
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FIG. 5: Decay of the energy spectra E{k) in the dissipation range for different Reynolds numbers, (a) Local slope of the 
spectrum as a function of the wavenumber k = krjK- The symbols for the different Taylor microscale Reynolds numbers are 
indicated in the legend of the figure. The spectra have been shifted with respect to each other for better visibility, (b) Local 
exponent aioc as a function of the wavenumber k. Least square fits to (|12p are done in the interval [k — 1/2, k + 1/2]. (c) Local 
exponent Pioc as a function of k. 



VI. CONCLUDING REMARKS 



The finest-scale intermittent fluctuations of fluid turbulence have been studied with a spectral resolution never 
appUed before. They are associated with a whole range of local dissipation scales rather than a mean dissipation 
scale - the Kolmogorov length rjK- We find that the increase of small-scale intermittency with increasing Reynolds 
number goes in line with an increasing spread-out of the local dissipation scales into the dissipation range. The 
growth of this scale range with respect to the Reynolds number is small, but present. This is consisitent with a 
logarithmic dependence of the extension of the intermediate dissipation range on the Reynolds number as proposed in 
01 . Furthermore, we detected a growing amplitude of sub-Kolmogorov-scale fluctuations which manifests in a slower 
exponential decay of the energy spectra in the large-wavenumber range. Both numerical results conflrm indications 
from other studies that increasing small-scale intermittency affects a growing number of scales in the dissipation range 
that are expected to be Reynolds-number-independent (and thus universal) in the classical theory of turbulence. 
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